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, ABSTRACT 
l> . 

, We construct hydrogen atmosphere models for magnetized neutron stars in radiative equilibrium with surface 
\ fields B = 10^^ — 5 X 10^^ G and efi'ective temperatures Teg ~ a few x 10^ — 10^ K by solving the full radiative 
transfer equations for both polarization modes in the magnetized hydrogen plasma. The atmospheres directly 
C\J determine the characteristics of thermal emission from isolated neutron stars. We study the effects of vacuum 
polarization and bound atoms on the atmosphere structure and spectra. For the lower magnetic field models 
{B ~ lO^'^ G), the spectral features due to neutral atoms lie at extreme UV and very soft X-ray energies and 
therefore are not likely to be observed. However, the continuum flux is also different from the fully ionized 
O ' case, especially at lower energies. For the higher magnetic field models, we find that vacuum polarization 
^ ■ softens the high energy tail of the thermal spectrum. We show that this depression of continuum fiux 
' strongly suppresses not only the proton cyclotron line but also spectral features due to bound species; 
therefore spectral lines or features in thermal radiation are more difficult to observe when the neutron star 
magnetic field is ^ 10^"^ G. 
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INTRODUCTION 

Thermal radiation from the surface of isolated neutron stars (NSs) can provide invaluable information on 
the physical properties and evolution of NSs. Such radiation has been detected in various types of isolated 
NSs: from radio pulsars (see Becker 2000; Pavlov et al. 2002) to old and young radio-quiet NSs (see Treves et 
al. 2000; Pavlov et al. 2002) to soft gamma-ray repeaters (SGRs) and anomalous X-ray pulsars (AXPs) (see 
Hurley 2000; Israel, Mereghetti, & Stella 2001; Mereghetti et al. 2002), which form a potentially new class of 
NSs ("magnetars") endowed with superstrong [B ^ 10^^ G) magnetic fields (see Thompson &: Duncan 1996; 
Thompson 2001). The NS surface emission is mediated by the thin atmospheric layer (with scale height 
~ 0.1 — 10 cm and density ~ 0.1 — 100 g/cva?) that covers the stellar surface. Therefore, to properly interpret 
the observations of NS surface emission and to provide accurate constraints on the physical properties of 
NSs, it is important to understand in detail the radiative properties of NS atmospheres in the presence of 
strong magnetic fields. 

Steady progress has been made over the years in modeling NS atmospheres (see Pavlov et al. 1995; Ho & 
Lai 2001, 2003; Zavlin & Pavlov 2002 for more detailed references on observations and on previous works of 
NS atmosphere modeling). These atmosphere models have played a valuable role in assessing the observed 
spectra of radio pulsars and radio-quiet NSs. We have studied the H and He atmospheres of NSs with 
magnetic fields B ~ 10^^ - 10^^ G and effective temperatures %s ~ 10^ - 10^ K (Ho & Lai 2001, 2003; 
Ho et al. 2002). Here we briefly discuss our atmosphere model and the effects of vacuum polarization and 
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report on recent work in incorporating hydrogen bound species in the models to produce partially ionized 
atmospheres. 

ATMOSPHERE MODEL 

We consider an isolated NS with a plane-parallel atmosphere. This is justified since the atmospheric 
scale height H ^ 10 cm is much less than the NS radius i? 10 km. The atmosphere is composed of pure 

hydrogen. A uniform magnetic field B, aligned perpendicular to the surface, permeates the atmosphere, and 
the spectra presented indicate emission from a local patch of the NS surface with the indicated magnetic 
field strength and effective temperature. 

In the highly magnetized plasma that characterizes NS atmospheres, there are two photon polarization 
modes, and they have very different radiative opacities. The two photon modes are the extraordinary mode 
(X-modc), which is mostly polarized perpendicular to the k — B plane, and the ordinary mode (0-mode), 
which is mostly polarized parallel to the k — B plane, where k specifies the direction of photon propagation 
and B is the direction of the external magnetic field. The 0-mode has a significant component of its electric 
field E along B for most directions of propagation, except when k is nearly parallel to B, and therefore the 
0-mode opacity is close to the B = value, while the X-mode opacity is much smaller. 

To construct self-consistent atmosphere models requires successive iterations, where the temperature pro- 
file T{t) is adjusted from the previous iteration in order to satisfy radiative equilibrium (r is the Thomson 
depth within the atmosphere). We solve the radiative transfer equations (RTEs) for the two coupled pho- 
ton polarization modes, together with the boundary conditions, by the finite difference scheme described 
in Mihalas (1978). We use a variation of the Unsold-Lucy temperature correction method as described 
in Mihalas (1978) but modified to account for full radiative transfer by two propagation modes in a mag- 
netic mcdivmi. The process of determining the radiation intensity from the RTE for a given temperature 
profile, estimating and applying the temperature correction, and then recalculating the radiation intensity 
is repeated until convergence of the solution is achieved (see Ho & Lai 2001 for details of our numerical 
method) . 

VACUUM POLARIZATION 

It is well-known that polarization of the vacuum due to virtual e'^e~ pairs becomes significant when 

B Bq, where Bq = m^(?/eh = 4.414 x 10^'^ G. Vacuum polarization modifies the dielectric prop- 
erty of the medium and the polarization of photon modes, thereby altering the radiative scattering and 
absorption opacities (see Ho & Lai 2003 and references therein). Of particular interest is the "vacuum 
resonance" phenomenon, which occurs when the effects of the vacuum and plasma on the linear polariza- 
tion of the modes cancel each other, giving rise to a resonance feature in the opacity at photon energy 
Ey ~ 1.02 (Z/yl)(/9/l g cm" 

-3)1/2(5/1014 G)-i/(5) keV, where Z and A are the atomic charge and mass 
of the ion, respectively, p is the density, and f{B) is a slowly-varying function of B of order unity. 

It was shown in Lai & Ho (2002, 2003) (see also Gnedin et al. 1978; Pavlov & Gnedin 1984) that high energy 
photons propagating in the atmospheric plasma can adiabatically convert from one polarization mode into 
another at the vacuum resonance density pv ^ 0.96 (A/Z)(£'/l kcV)^(i?/10"'^^ G)'^f{B)~'^ g cm"^. Across 
the resonance, the orientation of the polarization ellipse rotates by 90° , although the helicity does not change 
(see Figure 1). This resonant mode conversion is analogous to the Mikheyev-Smirnov-Wolfenstein (MSW) 
effect for neutrino oscillations. Because the two photon modes have vastly different opacities, this vacuum- 
induced mode conversion can significantly affect radiative transfer in strongly magnetized atmospheres. 
Figure 2 illustrates the main effect that vacuum polarization has on the atmosphere spectrum. The left side 
of Figure 2 shows the decoupling densities of the 0-mode and X-mode photons (i.e., the densities of their 
respective photospheres) when the vacuum polarization effect is neglected. The right side of Figure 2 shows 
that the effective decoupling depths of the photons are changed when the vacuum resonance lies between 
these two photospheres. We see from Figure 2 that mode conversion makes the effective decoupling density 
of X-mode photons (which carry the bulk of the thermal energy) smaller, thereby depleting the high energy 
tail of the spectrum and making the spectrum closer to blackbody (although the spectrum is still harder 
than blackbody because of non-grey opacities). This expectation is borne out in self-consistent atmosphere 



3 



lowp 

::<:::> 

0-mode 



V 



highp 



X-mode 



A 

\/ 

\j 

X-mode 







-> 



0-mode 



Fig. 1. A schematic diagram illustrating mode 
conversion due to vacuum polarization: As an 
X-mode (0-mode) photon at high density tra- 
verses the vacuum resonance density pv, it be- 
comes an 0-mode (X-mode) photon at low 
density if the adiabatic condition is satisfied. In 
this evolution, the polarization ellipse rotates 
90°, and the photon opacity changes signifi- 
cantly. 
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Fig. 2. A schematic diagram illustrating how 
vacuum polarization-induced mode conversion 
affects radiative transfer in a magnetar atmo- 
sphere. When the vacuum polarization ef- 
fect is turned off, the X-mode photosphere 
(where optical depth ~ 1) lies deeper than 
the 0-mode. With the vacuum polarization ef- 
fect included, the X-mode effectively decouples 
(emerges) from the atmosphere at the vacuum 
resonance, which lies at a lower density than 
the (original) X-mode photosphere. 



modeling presented in Ho & Lai 2003 and illustrated in Figure 3|^ for the case of B = 5 X 10^^ G and 
Teff = 5 X 10^ K. 

Another important effect of vacuum polarization on the spectrum, first noted in Ho & Lai 2003 and 
illustrated in Figure 3, is the suppression of proton cyclotron lines (and other spectral lines; see below). The 
physical origin for such line suppression is related to the depletion of continuum flux; this depletion makes 
the decoupling depths inside and outside the line similar. 



PARTIALLY IONIZED ATMOSPHERES 

A major complication when incorporating bound species in the atmosphere models arises from the strong 
coupling between the center-of-mass motion of the atom and the internal atomic structure. Potekhin & 
Chabrier (2003) have constructed the first thermodynamically consistent equation of state and opacity 
models of partially ionized hydrogen plasmas in strong magnetic fields. From the work of Potekhin &; 
Chabrier (2003), we obtain tables for the equation of state and absorption opacities of hydrogen, including 
the bound states. We interpolate between the table values to obtain densities and absorption opacities 
corresponding to the numerical grid of our models.^ Thus we are able to obtain partially ionized hydrogen 

^ All the models presented here use the corrected free-free absorption, which properly accounts for electron-ion collisions 
(Potekhin & Chabrier 2003; see also Ho et al. 2002). 

^ Note that we calculate the polarization vectors of the medium assuming the medium is fully ionized since the abundance of 
bound species is very low for the cases considered; this is not strictly correct (see Bulik & Pavlov 1996 for the case where the 
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atmosphere models (Ho et al. 2002). 

Figure 4 shows the spectra of hydrogen atmospheres with B = 10^^ G and T^ff = 5 x 10^ K, along 
with the non-magnetic hydrogen atmosphere model at the same Tefj and the blackbody spectrum with 
T = 5 X 10^ K. The proton cyclotron line at 6.3 eV is clear in both the fully ionized and partially ionized 
spectra.0 It is clear from Figure 4 that spectral features due to bound species at this field strength lie 
within the extreme UV (E ^ 50 eV) to very soft X-ray (E ^ 0.2 keV) regime and therefore are likely to 
be unobservable due to interstellar absorption. However, the effect of bound species on the temperature 
profile of the atmosphere and the continuum flux is significant. In particular, the optical flux is higher for 
the partially ionized atmospheres compared to the fully ionized atmospheres for the same Tefj. On the other 
hand, the partially and fully ionized models both yield very similar (neglecting the spectral features) X-ray 
flux for a given effective temperature, with the partially ionized model fluxes slightly lower than those from 
the fully ionized models; thus fitting only the high energy flux with fully ionized models would yield fairly 
accurate NS temperatures. 

Figure 5 compares the partially and fully ionized models (neglecting vacuum polarization) of hydrogen 
atmospheres with B = 5 x 10^^ G and T^g = 5 x 10^ K. Neglecting vacuum polarization allows us to 
distinguish features due to neutral atoms in the atmosphere; broad absorption features due to bound-free 
and bound-bound transitions at E ^ 0.76 keV and 4 keV, respectively, are apparent (the latter is from 
blending of two bound-bound transitions at 3.4 keV and 6.5 keV). The continuum flux between the models 
is not significantly different, and as noted in Ho & Lai (2001, 2003) and Zane et al. (2001), the proton 
cyclotron line is very broad when vacuum polarization is neglected. Figure 6 compares the models with 
vacuum polarization (complete mode conversion; see Ho &: Lai 2003). We see that vacuum polarization 
softens the high energy tail of the spectrum in the partially ionized model, like the fully ionized model. 
In addition, we see that vacuum polarization also suppresses spectral features due to the bound species 
of hydrogen. The reduced width of the proton cyclotron line and spectral features associated with bound 
transitions makes these features difficult to observe with current X-ray detectors, and we suggest that the 
absence of lines in the observed spectra of several AXPs (Patel et al. 2001; Juett et al. 2002; Tiengo et 
al. 2002) may be an indication of the vacuum polarization effect at work in these systems. 
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vacuum polarization. The solid line is for a par- 
tially ionized atmosphere, the dashed line is for 
a fully ionized atmosphere, and the dotted line 
is for a blackbody with T = 5 x 10^ K. 




E (keV) 



Fig. 6. Spectra of partially ionized hydro- 
gen atmospheres with B = 5 x 10^^ G and 
^eflf = 5 X 10^ K. The solid line is for an atmo- 
sphere which includes vacuum polarization, the 
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